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ABSTRACT 
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. Context. 

Aims. Solid O2 has been proposed as a possible reservoir for oxygen in dense clouds through freeze-out processes. The 
t-H ' aim of this work is to characterize quantitatively the physical processes that are involved in the desorption kinetics of 

CO-O2 ices by interpreting laboratory temperature programmed desorption (TPD) data. This information is used to 
t-H ■ simulate the behavior of CO-O2 ices under astrophysical conditions. 

■ Methods. The TPD spectra have been recorded under ultra high vacuum conditions for pure, layered and mixed mor- 

£N| ' phologies for different thicknesses, temperatures and mixing ratios. An empirical kinetic model is used to interpret the 

£NJ , results and to provide input parameters for astrophysical models. 

C<") ■ Results. Binding energies are determined for different ice morphologies. Independent of the ice morphology, the des- 

£SJ ' orption of O2 is found to follow th -order kinetics. Binding energies and temperature-dependent sticking probabilities 

for CO-CO, O2-O2 and CO-O2 are determined. O2 is slightly less volatile than CO, with binding energies of 912±15 
£-■■«. , versus 858±15 K for pure ices. In mixed and layered ices, CO does not co-desorb with O2 but its binding energies 

are slightly increased compared with pure ice whereas those for O2 are slightly decreased. Lower limits to the sticking 
*~^« ■ probabilities of CO and O2 are 0.9 and 0.85, respectively, at temperatures below 20 K. The balance between accretion 

' and desorption is studied for O2 and CO in astrophysically relevant scenarios. Only minor differences are found between 

SpH, the two species, i.e., both desorb between 16 and 18 K in typical environments around young stars. Thus, clouds with 

q . significant abundances of gaseous CO are unlikely to have large amounts of solid O2. 

S_i ' Conclusions. 
& • 
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1. Introduction dance in cold dark clouds in the range of 3x 10~ 6 -lx 10~ 7 

(Goldsmith et al. 2000, Pagani et al. 2003, Liseau et al. 
C$ • Gas-grain interactions play a key role in the chemical evolu- 2005). This low abundance, along with the low abundance 
' tion of star- forming regions. During the first stages of star of gaseous H 2 0, raises serious questions about the to- 
formation virtually all species accrete onto grains in dense tal oxygen budget when compared with the well observed 
cold cores. Later on in the star formation sequence — when atomic oxygen abundance of 3xl0~ 4 in diffuse clouds 
so-called hot cores are formed — grains are warmed to tern- (Meyer et al. 1998, Ehrenfreund & van Dishoeck 1998). 
peratures where molecules can desorb again. In order to One possible explanation for the 'missing' oxygen is that 
characterize this astrophysical process quantitatively it is it is frozen out onto grains in the coldest regions. The fun- 
necessary to understand the underlying molecular physics damental vibration of solid 2 around 1550 cm -1 (6.45 fim) 
by studying interstellar ice analogs under laboratory con- becomes observable through perturbations of the symme- 
trolled conditions. A series of recent papers shows that even try of 2 in a matrix or ice containing other molecules 
for simple molecules such a quantification is far from triv- (Ehrenfreund et al. 1992). This band has been sought to- 
ial (e.g. Fraser et al. 2001, Collings et al. 2003, 2004, Oberg wards the proto-stellar sources RCrA IRS2 and NGC 7538 
et al. 2005, Bisschop et al. 2006). In the present work re- IRS9 (Vandenbussche et al. 1999), but not detected. Upper 
suits for CO and 2 ices are discussed that extend recent limits between 50% and 100% of solid 2 relative to solid 
work comparing CO-N 2 and CO-0 2 ice features (Fuchs et CO have been reported from analysis of the CO profile, 
al. 2006) to an empirical kinetic model characterizing the since solid CO, in contrast with solid 2 , has been ob- 
desorption behavior. served through its vibrational band at 4.67 /im (2140 cm^ 1 ) 
The reason for focusing on ices containing 2 is that a ( e -§- Chiar et aL 1998 > Pontoppidan et al. 2003). Here the 
substantial amount of interstellar oxygen may well freeze amount of frozen 2 can be estimated by observing its mflu- 
out onto grains in the form of molecular oxygen. Attempts ence on the sha P e of the C0 absorption band. Transmission 
to determine gaseous 2 -abundances from recent SWAS s P ectra recorded for mixed CO-0 2 ices show indeed signif- 
ied ODIN campaigns put upper limits on the 2 abun- icant changes compared to pure CO ices (Ehrenfreund et 
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Fig. 1. The cryo-pumping factor (C X (T) = p r t/pt) as a 
function of temperature for CO (solid) and O2 (dashed). 

al. 1997, Elsila et al. 1997). However, such changes can also 
be caused by grain size and shape effects (Dartois 2006). 
The best limits on solid O2 therefore come from analysis 
of the weak solid 13 CO band which is not affected by grain 
shape effects. This band leads to upper limits of 100% on 
the amount of O2 that can be mixed with CO (Boogert et 
al. 2003, Pontoppidan et al. 2003). 

In this study, laboratory results for CO-O2 ices are pre- 
sented, both in pure, layered and mixed ice morphologies for 
varying ice thicknesses and for different relative abundances 
(Section 2). The focus is on temperature programmed des- 
orption (TPD) data that are recorded to study the des- 
orption process and to visualize changes in ice morphology 
during heating (Section 3). An empirical kinetic model is 
used to interpret the TPD data (Section 4). The aim of this 
work is to derive accurate molecular parameters (CO-CO, 
CO-O2 and O2-O2 binding energies, classification of des- 
orption kinetics, and temperature dependent sticking co- 
efficients) to allow reliable predictions of typical behavior 
under astrophysical conditions (Section 5). 

2. Experiment 

The experimental setup consists of an ultra-high vacuum 
setup (background pressure < 5xl0~ n mbar) in which in- 
terstellar ice analogs are grown with mono-layer precision 
on a 2.5x2.5 cm 2 sample surface made out of a 0.1 /im 
thick polycrystalline gold film. Physical-chemical interac- 
tions can be studied using TPD and/or reflection absorp- 
tion infrared spectroscopy (RAIRS). Details of the setup 
are available from van Brockhuizcn (2005) and Fuchs et al. 
(2006) and here only relevant details are given. We assume 
that for a sticking coefficient of 1, molecules hitting the cold 
surface (14 K) build up one monolayer (ML) at a steady 
gas exposure of 1.33-10 -6 mbar s _1 (10 -6 Torr s _1 ). This is 
by definition one Langmuir [L] and corresponds to roughly 
10 15 particles cm~ 2 in a non-porous condensed solid. In the 
present experiment amorphous ice is grown using a 0.01 L 
s _1 growth velocity resulting in a well specified layer thick- 



Table 1. Overview of all ice morphologies used in the CO- 
O2 experiments. The thickness is given in Langmuir [L] . 
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ness as long as the growth is linear with the exposure time. 
Typical layer thicknesses are generated in the range be- 
tween 20 and 80 L as astrophysical observations of young 
stellar objects in nearby low-mass star-forming clouds indi- 
cate that CO ices exist at thicknesses around 40 monolayers 
(Pontoppidan et al. 2003). In addition, astrochemical mod- 
els suggest that O2 may freeze out onto a pre-existing CO 
layer and consequently both pure, mixed and layered ices 
are studied here (Hasegawa et al. 1992, d'Hendecourt et al. 
1985, Bcrgin et al. 2000). Studies of 2 in an H 2 0-rich en- 
vironment have been performed by Collings et al. (2004). 

Two sets of ices have been investigated, 13 CO- 18 02 and 
12 CO- 16 02, both using isotopes of at least 99% purity. Most 
experiments were performed using 13 CO- 18 2 to distin- 
guish between possible impurities in the vacuum system. 
An overview of the used ice samples is given in Table 1. 
Throughout this paper, the notation x/y and x:y denote 
layered (x on top of y) and mixed ice morphologies, re- 
spectively. The notation 1/1 indicates layered ices of equal 
thickness, 1:1 refers to equally mixed ices and x/40 L stands 
for ices where the thickness of the top layer is varied by x L 
and the bottom layer is kept constant at 40 L. 

Once an ice is grown the desorption behavior is exam- 
ined by a controlled linear temperature rise from 12 to 80 
K. The heating rate in all experiments is 0.1 K min -1 , un- 
less stated differently. The CO and 2 molecules that are 
released are monitored by a quadrupole mass spectrome- 
ter (QMS). For a constant pumping speed the QMS sig- 
nal strengths are proportional to the amount of desorbed 
species. Since the deposition surface, together with other 
cold surface areas like the radiation shield, can be regarded 
as a cryogenic pump a temperature dependent cryo-pump 
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Fig. 2. TPD spectra of pure ices, (a) 13 CO and (b) 18 2 
- for exposures of 20 L (solid), 40 L (dot-dashed), 60 L 
(dashed) and 80 L (dotted) - and (c) 12 CO and (d) 16 2 - 
for exposures of 20 L (solid), 40 L (dot-dashed) and 80 L 
(dotted) . The grey vertical lines indicate the peak temper- 
atures at 40 L. 

factor has to be introduced as the pressure in the chamber 
decreases when the temperature drops below 50 K. This 
effect is measured by opening the leakage valve up to a 
certain pressure (p r t) at room temperature after which the 
system is cooled down to 14 K and the corresponding pres- 
sure (jpt) is determined. In Fig. 1 the ratio C X (T) =p rt /p t 
is shown as a function of temperature for both CO and O2 
revealing that at 14 K the cryo-pumping factor increases 
by a factor of about 7 with respect to the pumping speed, 
St p , of the turbo pump alone. This effect is important and 
will be explicitly taken into account in the kinetic model. 

Reflection- Absorption Infrared (RAIR) spectra have 
been taken simultaneously with the TPD data but are not 
included here. Spectra for selected CO-O2 experiments were 
presented and discussed in Fuchs et al. (2006); they provide 
additional information on the processes occurring in the ice 
during heating. 

3. Experimental results 

TPD spectra of pure and mixed/layered CO-O2 ices are 
shown in Figs. 2 and 3, respectively. 

Pure ices. In the case of pure 13 CO ice (Fig. 2a) and pure 
12 CO ice (Fig. 2c) the desorption starts around 26.5 K and 
peaks between 28 and 29 K. All leading edges clearly coin- 
cide and peak temperatures shift towards higher tempera- 
tures for higher thicknesses. This behavior is typical for a 
0"*-order desorption process (see §4.1) and corresponds to 
a desorption rate that is independent of the ice thickness 
and that remains constant until no molecules are left on the 
surface (see also Collings et al. 2003, Bisschop et al. 2006). 
The binding energies for both isotopes are similar and con- 
sequently nearly the same peak temperature is found. For 
pure 18 02 ice (Fig. 2b) and pure 16 02 ice (Fig. 2d) a similar 
behavior is found as for CO: all desorption curves have the 
same leading edges and the peak temperature slightly shifts 
for increasing thicknesses (from 29.5 to 31.3 K). The des- 



Fig. 3. TPD spectra for mixed/layered ices, (a) O2 trace of 
(5, 10, 20, 30 and 40 L) 18 O 2 /40 L 13 CO - differential layer, 
(b) (20, 40 and 80 L) 18 2 (solid line) and 13 CO (dashed 
line) for 1/1 O2/CO system. (The shoulder observed in the 
40 L measurement is an experimental artifact absent in 
other experiments), (c) (10, 20, 40, 60, 80 L) 18 2 (solid 
line) and 13 CO (dashed line) for 1:1 mixed ices, (d) (20, 
40 and 80 L) 16 2 (solid line) and 12 CO (dashed line) for 
the 1/1 O2/CO system. The grey vertical lines indicate 
the peak temperatures for pure O2 (solid) and pure CO 
(dashed) for 40 L. 



orption peak of O2 is clearly shifted to higher temperatures 
by ~2 K compared with CO. A small shift of peak temper- 
atures of less than 0.5 K is observed when comparing 18 02 
and 16 2 . 

Layered ices. In Fig. [3^ TPD spectra are shown for dif- 
ferent layers of 18 2 on top of 40 L 13 CO and in Fig. [5b 
(Hi) for 18 2 ( 16 2 ) on top of 13 CO ( 12 CO) for different 
1/1 configurations (i.e. 20 L 18 2 /20 L 13 CO, 40 L 16 2 
/ 40 L 12 CO, etc.). In all cases the O2 desorption behaves 
as a O^-order process and is very similar to that observed 
for pure O2 ice. This can be expected as O2 is less volatile 
than CO, i.e., at the desorbing temperatures of O2 there is 
little CO left to influence the O2 desorption process. The 
only noticeable difference compared to pure ices is that the 
peak temperature shifts to a lower value for O2 and to a 
higher temperature for CO by about 0.5 K. This indicates 
that in the layered ice systems the binding energies change; 
for CO the binding energy increases whereas for O2 it de- 
creases. With exception of the 20 L 13 CO spectrum, all lay- 
ered CO traces reveal a O^-order process as can be seen in 
Fig. 3b and Fig. 3d for the 12 CO. Potentially, the presence 
of O2 on top of the CO ice can change the desorption pro- 
cess for CO but RAIR spectra of layered ices do not change 
much with respect to pure ices. The lack of co-desorption 
of CO with the O2 suggests that the molecular interaction 
between these species is weak. 

Mixed ices. In Fig. 3c the TPD spectra are shown for dif- 
ferent layer thicknesses of 1:1 18 02: 13 CO mixed ices. The 
O2 desorption follows again a 0*' l -order process, but the 
band as a whole is slightly broadened by 10-15% compared 
to pure O2. In the 13 CO spectra a shoulder around 29.6 K 
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Fig. 4. Model TPD spectra of pure ice (solid lines) com- 
pared with laboratory data (dashed lines) for 20, 40, 60 
and 80 L ices, (a) 13 CO, (b) 18 2 , (c) Model 13 CO (dashed) 
and 18 O2 (solid) using recommended energy values and (d) 
variation of x 2 with energy for an exposure of 40 L, illus- 
trating the well-determined minima. The effect of fits with 
and without using the cryo-pumping function C X (T) are 
shown. 



appears, i.e. at the O2 desorption temperature, which sug- 
gests that a fraction of the CO desorbs from CO-CO bind- 
ing sites like in pure ices whereas the rest desorbs from a 
mixed environment, e.g. CO-O2 binding sites. The main iso- 
topes (not shown in the figure) exhibit a similar behavior. 
Compared with previous experiments performed on CO-N2 
ices (Oberg et al. 2005, Bisschop et al. 2006) these results 
have been interpreted as follows (Fuchs et al. 2006): neither 
N 2 nor O2 possess an electric dipole moment so they inter- 
act with CO mainly via quadrupole interactions. However, 
solid O2 has a 4 to 6 times weaker quadrupole moment 
compared to N 2 and CO. Furthermore, N 2 and CO possess 
the same a-crystalline structure below 30 K, but a-phase 
2 has a different crystalline structure and also undergoes 
a phase change to the /3-form at 23.5 K. The combination 
of these two effects can lead to the absence of mixing and 
co-desorption in the CO-O2 system compared with CO-N 2 , 
as observed in our experiments. 



4. Empirical kinetic model of CO-O2 desorption 

The experimental TPD results are interpreted in terms of 
an empirical kinetic model, describing the desorption kinet- 
ics and providing values for fundamental molecular proper- 
ties to be used in astrochemical models. 



4.1. The model 

The kinetic desorption process for a species X can be ex- 
pressed by the well known Polanyi-Wigner type equations 
of the form, 



Rdes — 



dNg{X) 
dt 



= Vi[N s (X)} 1 exp[- 



Fig. 5. Model TPD spectra of layered and mixed ices; (a) 
Experimental (dashed) and model (solid) results of 13 CO 
in mixed ice; (b) Best fitting 13 CO and 18 2 model results 
for mixed ices; (c) Best fitting 18 02 model results for x/40 
L 18 2 / 13 CO and (d) best fitting model results for 13 CO 
and 18 O2 in 1/1 layered ices. 

with Rdes the desorption rate (molecules cm~ 2 s -1 ), N g (X) 
the number density (cm -2 ) of molecules desorbing from 
the substrate, Vi a pre-exponential factor (molecules 1-4 
cm 2 ( l_1 )s _1 ) f or desorption order i, N S (X) the number den- 
sity (cm -2 ) of molecules on the surface at a given time t, 
Ed{X) the binding energy (in K) and T the surface temper- 
ature (also in K). The desorption order reflects the nature 
of the desorption process and is expressed in integer values 
although also non-integer values arc allowed. 

Ideally, the pre-exponential factor and the desorption 
energy are determined independently. However, since both 
parameters have a similar effect on the fitting routine a 
unique parameter set could not be found. Therefore, the 
pre-exponential factor is approximated by the harmonic os- 
cillator of a solid in terms of a vibrational frequency (s -1 ) 

by 



2N s E d (X) 



(2) 



with N s 



10 



15 



7T 2 M 

cm -2 and M the mass of species X 
(Hasegawa et al. 1992). This equation gives values around 
10 12 s _1 . In previous papers (see e.g., Bisschop et al. 2006, 
Collings et al. 2003) the pre-exponential factor has been 
taken as a free parameter. In the present work this pa- 
rameter has been linked to Ed{X) following Eq.(2). For a 
tft, -order desorption process, v has been multiplied with 
the surface density N s , i.e. vq = v ■ N s = v ■ N s (t = 0). For 
a l st -order process the vibrational frequency is taken as a 
pre-exponential factor, i.e. v\ = v. Thus the only parameter 
that is floating in the present model fit to the experimental 
data is the binding energy Ed{X) of the species. 

Starting from Eq. (1), the kinetics are represented in terms 
of two coupled differential equations 



dN s (X) 
dt 



(1) 



J2 v t [N s (X)] l cx V [ 

i=0,l 



E d ,{X), 



(3) 
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Type of ice 


Reaction 


Rate equation 


V 


E d 


1 










[molecules 














cm a(i-1) s _1 ] 


[K] 


















Pure 


CO(s) 


- CO(g) 


^0 exp(-£ d /T) 


7.2E26 


858±15 





Pure 


2 (s) 


- o 2 (g) 


z/ exp(-£ d /T) 


6.9E26 


912±15 





Layered 


CO(s) 


-» CO(g) 


1/0 exp(-£ d /T) 


7.2E26 


856±15 





Layered 


2 (s) 


-» 02(g) 


i^o exp(-£ d /T) 


6.9E26 


904±15 





Mixed 


CO(s) 


- CO(g) 


z/i[CO] exp(-£ d /T) 


7.6E11 


955±18 


1 




CO(s) 


-> CO(g) 


1/0 exp(-£ d /T) 


7.2E26 


865±18 





Mixed 


2 (s) 


- 02(g) 


^0 exp{-E d /T) 


6.8E26 


896±18 



















Pure 


CO(s) 


- CO(g) 


i/ exp(-£ d /T) 


7.0E26 


854±10 





Pure 


2 (s) 


- o 2 (g) 


l/ exp(-£ d /T) 


6.6E26 


925±10 





Layered 


CO(s) 


- CO(g) 


1/0 exp(-£ d /T) 


7.1E26 


860±15* 





Layered 


2 (s) 


- 02(g) 


^0 exp(-£ d /T) 


6.5E26 


915±10 





Mixed 


CO(s) 


-» CO(g) 


i/i[CO] exp(-£ d /T) 


7.5E11 


965±10 


1 




CO(s) 


-f CO(g) 


i/ exp(-£ d /T) 


7.2E26 


890±10 





Mixed 


2 (s) 


- O a (g) 


^0 ex P (-£; d /r) 


6.2E26 


915±10 





Pump 
















CO(g) -> CO (pump) 


k pump [CO(g)] 


0.00024 




1 




2 (g) -> 2 (pump) 


kp Um p [02(g)] 


0.00036 




1 



* Average value for ice thicknesses greater than 30 L. For ices between 20 and 30 L use 888±15 K. 



»=0,1 

-C x {T)k„N g (X) (4) 

representing the density change of the solid phase (s) and 
gas phase (g) molecules. The last term of Eq. (Q} gives the 
removal of gaseous species from the vacuum chamber by 
the pump. Here, k pump is the pumping constant in s _1 and 
C X (T) is a dimensionless cryo-pumping factor between 1 
and 7 as discussed in Sj2l The pumping constant, k pump , is 
constrained by fitting the pump down curve of both species, 
i.e., the slope of the TPD curve at temperatures higher than 
the desorption peak temperature. In most cases, the peak 
tail is not well reproduced due to the presence of other cold 
surfaces in the system. Very small amounts of gas that have 
missed the target can get deposited and can be desorbed 
afterwards at a temperature different from that of the tar- 
get surface. A direct consequence is a deviation in model 
curves from the experimental plots at high temperatures. 
However, this does not affect the determination of the bind- 
ing energies that are mainly sensitive to the peak value and 
the leading edge. 

In order to calculate the temperature dependent rate 
measured in the TPD experiments, dn/dt is written as 

dn dn dT dn 

~dt ~ dT~dt = ^~dT ^' 

with dn/dT the temperature dependent rate (molecules 
cm _2 K~ 1 ), and (3 = dT/dt the constant TPD heating rate 
(K s _1 ) that corresponds to 0.1 K/minute in the present 
study unless stated differently 

To extract the binding energy from the observations a stan- 
dard minimization of \ 2 is used which represents the sum 



over the squares of the differences between the experimen- 
tal points and the calculated ones. In this procedure first 
a th - or l s *-order process is fitted and subsequently the 
appropriate parameters are optimized. 

4.2. Results 

The results are shown in Figs. [Hand [5] and Table 2 summa- 
rizes the model equations with best fit parameters. In nearly 
all experiments a linear dependence between ice thicknesses 
and desorption energies has been found. 

Pure ices. For pure 18 C>2 five experiments are performed 
using different thicknesses (20 L, 40 L, 50 L, 60 L, 80 L ) 
and the best fits to these experiments (Fig. UJd) give 932, 
927, 924, 928 and 918 K, respectively, corresponding to an 
average value of 925 ± 10 K. The error of the mean value is 
chosen to give a conservative estimate of these parameters. 
Similarly, for pure CO (Fig. |4^) the binding energy is 854 
± 10 K which is consistent with the value of 855 K reported 
by Oberg et al. (2005) and Bisschop et al. (2006). Both pure 
CO and pure O2 ice exhibit O'^-order kinetics. The TPD 
spectra arc reproduced very well with this model using only 
one free parameter as is demonstrated in Figs. [4ji and b, 
where the experimental and model results are plotted on 
top of each other. In Fig. 0Ji the variations of x 2 f° r the 
experiments of 40 L for the pure CO and O2 are shown. 
The energy values for which \ 2 l& minimum are taken as 
the final values. This yields 854 K for 13 CO and 927 K for 
18 02- As mentioned in the experimental section we have 
explicitly taken into account the effect of cryo-pumping. A 
systematic deviation of 2 to 3 K in the desorption energy is 
induced when this factor is neglected. This is illustrated in 
Fig. |IJi. Inclusion of the C X (T) function generally improves 
the fit substantially. 

12 CO and 16 02 show a similar desorption behavior and the 
binding energies are 858 ± 15 K and 912 ± 15 K, respec- 
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tively. Thus 12 CO has the same binding energy as 13 CO 
within the errors. For 16 02 ices the binding energy is low- 
ered by 1.6% with respect to 18 C>2 on average. Because of 
the smaller data set for the main isotope the errors in the 
fitted parameters of the 12 CO and 16 O2 isotopes are larger 
than those for 13 CO and 18 2 . 

Layered ices. The desorption of 18 02 in 1/1 (Fig. [5]i) 
and x/40 (Fig. [5};) layered ices is O t,l -order and the binding 
energy is determined as 915 ± 10 K, which is slightly lower 
than for the pure ices. 16 O2 shows a similar behavior with 
a binding energy of 904 ± 10 K. 13 CO desorption is a th - 
order process and only at low coverages a l s *-order process 
may be involved. Since there is no signature in the TPD 
spectra (nor in the RAIR spectra) of mixing or segregation 
and since the spectrum looks very similar to pure ices it is 
modeled in O'^-order with an average binding energy of 860 
± 15 K for ices thicker than 30 L. The correctness of this 
procedure is confirmed by the TPD spectra of 12 CO which 
exhibit a 0*' l -order desorption with a binding energy of 856 
± 15 K. 

Mixed ices. O2 in mixed ices (Fig. is again t/l -order 
with a similar binding energy as in layered ices. In mixed 
ices, TPD spectra of 13 CO are broader with respect to pure 
ices and a shoulder around the desorption temperature of 
O2 is observed. This suggests that CO is desorbing from 
a wider range of binding sites and this is also supported 
by RAIR spectra that are broader with respect to pure ice 
spectra. Both the peak and shoulder are fitted very well 
using a combination of th - and l st -order processes (see 
Fig. [Sh,). About 50% of the CO desorbs from CO-CO bind- 
ing sites as in the pure ices but with a binding energy of 
890 ± 10 K. The residual CO molecules desorb through a 
l st -order process from a mixed environment, i.e. including 
CO-O2 binding sites, with a binding energy of 965 ± 5 K. 

Table 2 has separate entries for the 13 CO- 18 02 and 
12 CO- 16 02 isotopomers. Independent of the used iso- 
topomers it can be concluded that the binding energies 
of CO increase in the following order E purc < -Eiaycrcd < 
-Emixcd- For 2 this is inverted E purc > E laycicd > E mix . 

The effect of (3 on the binding energies. In order to 
rule out any dependencies on the adopted heating rates, 
experiments have been performed for /3=0.1, 0.2 and 0.5 K 
min -1 on 40 L CO ices. This is of relevance when applying 
laboratory values to interstellar warm-up time scales. The 
calculated v and Ed values for the 0.5 K min -1 experiment 
have been used to predict the experimental TPD curve for 

0. 2 and 0.1 K min -1 . The deviation between the calculated 
and experimental values are within the experimental error, 

1. e. AEd = ± 15 K, with a slight tendency for lower inferred 
Ed at lower heating rates. 

4.3. Sticking probability 

In addition to the desorption rates, the sticking coefficients 
also play an important role to describe freeze-out onto in- 
terstellar grains. The measurement procedure has been ex- 
tensively discussed in Bisschop et al. (2006) and Fuchs et 
al. (2006). It is important to note that the measurements 
only provide an 'uptake' coefficient 7 rather than a stick- 
ing coefficient, S, as only the net rate of molecules sticking 
and leaving the surface can be given. Consequently, the 
values given in Fig. [5] represent lower limits for the sticking 




Temperature /K 



Fig. 6. Uptake coefficient 7 representing the lower limits 
on the sticking coefficients as functions of the temperature 
below 25 K. 

coefficients of O2 and CO. It is found that the freeze-out 
dominates for O2 on 2 , CO on CO and O2 on CO up to 
25 K with lower limits on the sticking probabilities between 
0.85 and 0.9, i.e., close to unity. Under real astrophysical 
conditions this value will increase and approaches unity at 
10 K. 



5. Astrophysical implications 

The scenarios put forward to explain the absence of gaseous 
O2 in interstellar clouds can roughly be divided into two 
categories: time-dependent models of cold cores invok- 
ing freeze-out of oxygen (in all its forms) onto grains 
(e.g., Bergin et al. 2000, Aikawa et al. 2005), and depth- 
dependent models of large-scale warm clouds invoking deep 
penetration by UV radiation in a clumpy structure enhanc- 
ing the O2 photodissociation (e.g., Spaans & van Dishoeck 
2002). The laboratory experiments presented here are rel- 
evant to the first scenario. In these models, atomic O is 
gradually transformed with time into O2 in the gas phase. 
Freeze-out of oxygen occurs on a similar timescale, and 
any O is assumed to be turned effectively into H 2 ice 
on the grains where it subsequently sticks. No formation 
of solid O2 on the grain is expected because of the pres- 
ence of atomic hydrogen, which is much more mobile at low 
temperatures. The H2O ice formation lowers the gaseous 
[0]/[C] ratio and potentially leaves only a small abundance 
of O2 in the gas. The amount of solid O2 and the remain- 
ing fraction of gaseous O2 of ~ 10 -8 — 10~ 7 with respect to 
H2 depends sensitively on the temperature at which O2 is 
frozen out. A related question is to what extent O2 differs 
in this respect from CO, since CO is readily observed in the 
gas phase and in solid form. 

Astrochemical models usually assume l st -order desorp- 
tion with rates in cm~ 3 s . Appendix A (see online mate- 
rial) summarizes the equations used to apply our laboratory 
results to astrochemical models for both I s * and th order 
kinetics. The latter is more appropriate for thick ices in in- 
terstellar clouds and consistent with our laboratory data. 
We consider the simple case of pure desorption and accre- 
tion of 12 CO and 16 02, without any gas-phase or solid-state 
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Fig. 7. Astrophysical simulations of gaseous CO and O2 abundances relative to the total gas + solid abundances. Panel 
(a) shows pure CO and O2 in steady-state, panel (b) 0* h -order desorption for heating rates of 1 K/1000 yr, panel (c) 
the same but for 1 K/10 6 yr, and panel (d) shows the desorption of CO from mixed (layered) ices for 1 K/10 3 yr. In 
the panels (b)-(d) the solid line indicates accretion corresponding to n(H2) = 10 3 cm -3 , the dashed line corresponds to 
accretion for n(fi2)=10 7 cm~ 3 . For the panels (b)-(d) the adopted initial gaseous abundances of CO is 10~ 4 and for O2 
it is 10 -4 and 10~ 6 with respect to H2. 



chemistry. Thus, the main difference with the equations in 
§4.1 is that desorption is now balanced by a freeze-out (ac- 
cretion) term, rather than a pumping factor. Also, the lab- 
oratory heating rate can be replaced by a heating rate ap- 
propriate for the astrophysical object under consideration, 
for example a protostar heating its envelope. 

Our model only considers classical dust grains with a 
radius of 0.1 /im. The density n(H2) is varied to simulate 
the effects for clouds of different densities. In th order, the 
results depend on the initial number densities of the CO and 
O2 molecules. For CO, an abundance 10 -4 relative to H2 is 
chosen, initially all on the grains. For O2, two options are 
considered. The first option is that most oxygen has been 
converted into O2 at an abundance of 10 -4 , as found in gas- 
phase models and the maximum allowed by the observed 
upper limits on solid 02- The second case considers a much 
lower O2 abundance of 10 -6 , as found in gas-grain models. 

Fig. EH shows the results for pure CO and pure O2 ices 
in the simplest case of steady-state. The temperature at 
which these species desorb depends strongly on density: 
the higher the density, the faster the accretion rate which 
needs to be balanced by desorption. Because of the larger 



binding energy of O2, it always desorbs at ~2 K higher 
temperatures than CO. 

Fig. Eb shows the result for th order desorption at a 
heating rate typical for a proto-stellar environment of 1 
K/1000 yr (Lee et al. 2005). For abundances of CO and 
O2 of 10 -4 , both species desorb at ~2.5 K higher temper- 
atures than the steady-state case at low densities; at high 
densities, there is little difference. If the O2 abundance is 
lowered to 10 -6 the th order curve shifts to lower temper- 
ature by ~2 K. Thus, in contrast with the steady-state or 
equal abundance cases, O2 can desorb at lower tempera- 
tures than CO if its abundance is significantly lower. This 
reversal is a specific feature of the 0*' 1 order desorption and 
is not found in the 1 st order formulation (see Appendix A in 
the online version) . An important question is whether such 
a situation is astrochemically relevant. For O2 abundances 
as low as 10 -6 , the coverage becomes less than a mono-layer 
and I s * order kinetics or other effects due to the peculiari- 
ties of the ice (e.g. polar vs. apolar environment, compact 
vs. porous ice) will determine the desorption behavior. 

Fig. Et shows the abundance curve of the species for 
a heating rate of 1 K/10 6 yr, as expected for a cold core 
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at near constant temperature. Since the time scale is in- 
creased by 3 orders of magnitude with respect to the previ- 
ously considered heating rate there is simply more time for 
the molecules to desorb and consequently the entire profile 
shifts to lower temperatures by 0.5 - 2.5 K depending on 
the density of the species. 

Finally, Fig. [Tji shows the abundance curves of CO and 
O2 from mixed ices. The graphs of mixed and layered ices 
are nearly identical, thus the layered ices are not shown sep- 
arately. The relative abundance curve is shifted to a higher 
temperature by about 0.5 K compared with the pure ices 
because of the slightly higher binding energies. 

Overall, Fig. [7] shows that the differences in the des- 
orption behavior of O2 and CO with temperature are very 
minor for a wide range of realistic cloud densities and abun- 
dances. Thus, it is unlikely that a large reservoir of solid 
O2 is hidden in the bulk of molecular clouds which show 
abundant gaseous CO but no O2, unless the O2 is in a 
more strongly bound ice environment. Conversely, any re- 
gion with significant CO freeze-out should also have some 
solid O2. As noted in the introduction, the best limits come 
from analysis of the weak solid 13 CO band, which gives 
upper limits of 100% on the amount of O2 that can be 
mixed with CO, i.e., about (0.5 — 1) X 10~ 4 with respect to 
H2 (Boogert et al. 2003, Pontoppidan et al. 2003). Direct 
freeze-out of the gas-phase O2 abundances inferred from 
the ODIN measurements would give much lower limits. 

The small differences between CO and O2 desorption 
found here may become relevant in the interpretation of 
high spatial resolution observations of individual cold cores 
with temperatures in the 10-20 K range such as the pre- 
stellar core B68. 

6. Concluding remarks 

The desorption processes of CO-O2 pure, mixed and lay- 
ered ice systems have been investigated experimentally and 
modeled using an empirical kinetic model. The resulting 
molecular parameters can be used to model the desorption 
behavior of these ices under astrophysical conditions. We 
find that both pure 16 02 and pure 12 CO desorb through 
O^-order processes with binding energies of 912 ± 15 K 
and 858 ± 15 K, respectively. In mixed and layered ices the 
16 02 binding energy decreases to a lower value around 896 
± 18 K and 904 ± 15 K respectively. The 12 CO desorption 
from layered ices is 0*^-order with a binding energy of 856 
± 15 K. In mixed ices a combination of th - and l s *-order 
is found with desorption energies of 865 ± 18 K and 955 
± 18 K, respectively. For 18 02 and 13 CO, these numbers 
change by a few percent. 

O2 is less volatile than CO but CO does not co-desorb 
with O2. This is in contrast with the CO-N2 ice system for 
which Bisschop et al. (2006) found that N2 is more volatile 
than CO and that significant amounts of N2 co-desorb with 
CO. The sticking coefficients of CO and O2 at temperatures 
below 20 K are close to unity, with 0.85 as a lower limit. 
In cold clouds (T4 < 18 K), O2 can be frozen out onto the 
grains, but the relative difference in desorption between CO 
and O2 is so small that this is unlikely to be the explanation 
for the missing gaseous O2 in interstellar clouds which show 
significant gaseous CO. 
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Appendix A: Modeling accretion and desorption in 
astrochemical applications 

The gas-phase density n g (X) (in cm~ 3 ) of species X with 
respect to its total gas + solid density n to t(X) — n g (X) + 
n s (X) can be calculated from the accretion and desorption 
rates. For reasons of simplicity and to compare with other 
work we consider here the simplest case of a single type of 
grain with a classical radius r c —0.1 /im, 10 15 sites per cm 2 
and a grain abundance n sr /n(H2) = lCP 12 . 



A.l. Accretion rate 

The most general formulation of the accretion rate in cm 
s _1 is 

R acc = aSvn g {X)n gr (A.l) 

where er=7rr 2 is the grain cross section, S is the sticking 
coefficient (taken to be unity at T gr < 20 K) and v is the 
mean speed of the gas in cm s . For a classical grain and 
the canonical grain abundance this becomes 



R a 



4.55 x 10 



' i iwx)T ns(x)nH ' <a - 2) 

where T is the gas temperature and M(X) is the mass 
of X in atomic mass units (amu). The accretion rate per 
molecule X in s _1 is denoted as 

X = Racc/n g (X) (A.3) 

and the timescale for freeze-out can be computed from 
r = 1/A. Note that in dense cold cores, the grains may 
have grown to larger /im sizes than assumed here due to 
coagulation, in which case the abundance has to be low- 
ered. Also, a more sophisticated calculation should include 
a grain size distribution. 

A. 2. Desorption rate 

The desorption rate in the laboratory surface science exper- 
iments is given by Eq. (1) in cm -2 s , which is rewritten 
here as 



Rde 



dab 



Vi N % exp 



Jk 
kT, 



yr 



(A.4) 



where TV has the units of molecules cm 2 , and the unit of 
Vi changes with the order i of desorption as molecules 1-1 



cm 



2(i-l) 



For -order vq is in molecules cm s 



and for l st -order v\ is in s 1 . Ei is the binding energy for 
order i. 



A.2.1. I s * order 

For I s * order kinetics, the translation from the laboratory 
data to a astronomical desorption rate £ in s _1 is straight- 
forward: 

£1 = vi exp 



kT n , 



(A.5) 



with v\ and E\ taken directly from the laboratory data. 
The rate Rdes in enr 3 s _1 becomes 

Rdes = £in s {X) = iy 1 n s (X)exp 



(A.6) 



with n s (in cm 3 ) the density of molecules on the grain 
surface. 



A.2.2. 0*' 1 order 

For th order kinetics, the desorption rate Rdes in cm~ 3 s _1 
should be proportional to the number density of grains n gr 
and the average surface area A gr of a grain. The rate 'per 
grain' in s _1 



£0 = voA gr exp(-E a /kT gr ) 



(A.7) 

with i/q and Eq derived from the laboratory data for th 
order. The total rate is 



Rde 



£,on gr = u A gr n gr exp(-E /kT gr ) 



(A.8) 



A.3. Steady-state model 



If no chemical reactions in the gas or on the grains are 
included, a steady-state will be reached where the accretion 
balances the desorption. In 1 st order, this becomes 

Xn g (X) = Sm a (X) (A.9) 

with £1 from Eq. IA.5I Thus the ratio of gas-phase to solid- 
state molecules is given by 



^ = ii=2.2. 
n., A 



10 17 ^ 



'M 



exp(-Ei/kT gr ) 



T 



(A.10) 



Results of Eq lA.101 arc shown in Fig. lA.ll for two densi- 
ties n(H2) = 10 3 and 10 7 cm~ 3 . No heating rate j3 is involved 
in this calculation. It is the most simple "astrochemical" 
modeling result achievable. One note of caution, however: 
it is not correct to use laboratory values for Ei which were 
fitted to O'^-order laboratory data for l s *-order astrochem- 
ical equations and just change v. If I s * order kinetics are 
used (even though the laboratory data indicate th order), 
a better approach would be to first fit the laboratory data 
with 1 st order kinetics and then apply those in the astro- 
chemical models. The differences in binding energies can be 
as large as 12% for CO and O2. For pure N 2 , CO and O2, 
1 st order fits to the peaks of our laboratory curves give 800 
K, 955 K and 1035 K, respectively, with v x = 1.0 • 10 11 s" 1 
for N 2 (Bisschop et al. 2006) and v x = 1.0- 10 12 s" 1 for CO 
and O2. Figures [Th. and IA. II use these values. 

A.4. Time-dependent model 
A.4.1. 1 st order 

In the I s * order approach, the desorption rate changes with 
the number of species X bound to the grain, which is jus- 
tified only for (sub)monolayer or irregular coverage. In this 
scenario the abundance of a species is determined by two 
differential equations and some boundary (initial) condi- 
tions 



dn g /dt 
dn s /dt - 



-- -\n g + Cin s 
+Xn g — £xn s . 



(A.ll) 



These are the equations commonly used in time- 
dependent astrochemical models (e.g., Bergin & Langer 
1997, Charnley et al. 2001, Aikawa et al. 2005). 

For protostellar cores in which a star turns on, the sim- 
plest assumption is that there is a linear relation between 
the temperature T gr and the time t 



T gr (t) — T gr fi + dT gr /dt ■ t — T gr fi + fit 



(A.12) 
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with T gr! o around 10 K. At t = 0, all molecules are as- 
sumed to be frozen out onto grains, n s (t = 0) = n to t(X), 
as appropriate for a pre-stellar core prior to star formation. 
A completely opposite situation, not studied here, would 
be a quiescent cloud where dT gr /dt is close to zero and a 
molecules are initially in the gas phase. Note that in th 
case the timescales are dominated by the freeze-out proces 
rather than the evaporation timescale, and that this resull 
in a very different time behavior, depending on tempen 
ture. 

Fig. IA.1I includes the time-dependent I s * order resull 
for a heating rate of 1 K per 10 3 yr and two densitie: 
Compared with the steady-state results, the low densit 
curves are shifted to higher temperatures by ~3 K, whereE 
the high density curves are nearly identical. A slower heai 
ing rate of 1 K per 10 6 yr brings the low density curve 
closer to the low density steady-state results. 



A.4.2. th order 

In the th order approach, the desorption rate is constar 
with time, independent of the surface density n s (X). Th 
time-dependent equations now become 



-\n g 
-Xn g 



dng/dt — 
dn s /dt = +\n g — Rdes for n s > 



Rdes for n s > 
for n s = 



(A.l; 



= +\n a 



for n s = 



with Rdes from Eq. IA.8I This formulation is most appropr 
ate for thick, fully covered ice layers. For partial coveragi 
Rdes would need to be multiplied by a relative surface occi 
pation number 6 between and 1, but for the cases of pui 
CO and O2 ices considered here 8=1 is taken. The cond 
tion on the right hand side causes an abrupt end to th 
desorption process once the molecules on the surface have 
been evaporated. Instead of getting a smooth S— shaped 
curve like in the l st st-order approach, the density curve of 
a O'^-order formulation reveals a sharp edge at the end of 
the evaporation process. 

The results now depend on the choice of the initial 
n s (X)/n(H 2 ) abundance ratio of species X. For CO, a log- 
ical choice is n g (CO)/n(H2) = 10~ 4 . For O2, an extreme 
case would be to put all oxygen into O2, n ff (02)/n(H2) > 
10~ 4 . However, models and observations (see §1) sug- 
gest that a more plausible abundance is at least two or- 
ders of magnitude lower. Figure I A. 1 1 shows the results for 
n g (0 2 )/n(H 2 ) = 1(T 6 and 10~ 4 . The th order curves for 
CO and O2 for 10~ 4 differ by less than 0.7 K from those 
at I s * order. However, those for O2 at 10~ 6 are shifted to 
lower temperature by ~2.5 K. 

These examples illustrate the importance of proper 
modeling of both the order of the desorption processes in- 
volved and the balance with accretion. Shifts of a few K 
compared to the steady-state solution can occur either way 
depending on cloud parameters. Such small shifts are un- 
likely to affect the interpretation of large-scale molecular 
clouds, but they may potentially become significant in in- 
terpretations of the relative behavior of CO, O2 and N2 in 
cold dense clouds in the 10-20 K range. 

Finally, it should be noted that our results using the em- 
pirical model for pure ices are consistent with first principle 



considerations using saturated vapor pressure equations, as 
done for pure CO by Leger (1983). 
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Fig. A.l. The gas-phase fraction of CO (top) and O2 (bot- 
tom) for densities of n(H 2 ) = 10 3 cm~ 3 (full lines) and 10 7 
cm -3 (dashed lines) in the case of pure ices. Both steady- 
state (black), O t?l -order (blue/green) and l st -order (red) 
desorption are considered. In th order, the adopted abun- 
dance for CO is 10~ 4 whereas for O2 both 10~ 4 (green) and 
10~ 6 (blue) are considered. 



